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ABSTRACT 

Using ultra-precise data from space instrumentation we found that the underlying functions of stellar light curves from some AF pul¬ 
sating stars are non-analytic, and consequently their Fourier expansion is not guaranteed. This result demonstrates that periodograms 
do not provide a mathematically consistent estimator of the frequency content for this kind of variable stars. More importantly, this 
constitutes the first counterexample against the current paradigm which considers that any physical process is described by a contin¬ 
uous (band-limited) function that is infinitely differentiable. 

Key words. Asteroseismology - Methods: data analysis - Stars: oscillations - Kepler - 


1. Introduction 


The unexpe cted huge num ber of frequencies found in multi pe¬ 
riodic stars dPoretti |2009|) and the ubiquitous presence of corre¬ 
lated noise in the residuals of the fitting of the light curves of 
pulsating stars, could have a common origin. 

The necessary condition for a correct Fourier analysis, i.e. 
to have a Fourier expansion, is guaranteed when the function 
is analytic, otherwise Fourier analysis is not a consistent ap¬ 
proximation. In this work we examine the differentiability of the 
function des cribing photometric da ta of pulsatin g stars obtained 
by C oRoT (lAuverene et alJ l2009l) and Kepler loilliland et al.l 
20101) . We determine how smoothly the data points of the light 
curve can be htted in order to fully reproduce the function. We 
call this hne structure property the connectivity of the function 
at a given point. We studied this property by means of two nu¬ 
merical approac hes: cubic s plines - analyticQ, and autoregressive 
moving average (lBoxlll97^ methods - non-analytic. 


. Analyticity 


Fourier analysis is a useful technique for frequency detection. In 
order to apply Fourier te chniques to a giyen function, Parseyal’s 
theorem (iParseyalll 19921) must hold. This theorem states that the 
integral of the squared modulus of a function is equal to the inte¬ 
gral of the squared modulus of its Fourier transform. For this to 
be true the function must be square integrable. The Fourier ex¬ 
pansion of the function con yerges pointwise almost eyerywhere 
when it is square integrable (ICarlesonlll966l) . Unlike the Fourier 
series, which is an inhnite sum, the Discrete Fourier Transform 
(DFT) and its inyerse, the IDFT, which is the analogous to the 
Fourier series, always conyerge if the sequence is hnite, because 


* By definition only the first and second derivatives are continuous for 
cubic splines but we are interested here only in the first derivative. In 
any case, if the first derivative is not continuous the functions are non- 
analytic. 


they are finite sums. But when the Fourier series of the function 
is divergent, or not convergent to the value of the function at 
this point, the DFT no longer provides a mathematical descrip¬ 
tion based on the frequency components of the function. That is, 
a frequency is only a well-defined physical variable when the 
Fourier series converges. Therefore, before a Fourier analysis 
can be applied to a time series, the convergence of the Fourier 
expansion of the function must be demonstrated. 

On the other hand, the analyticity of a functio n is a sufficient 
condition for the square integrability (iDiik 11200^ . A function is 
said to be analytic if it is infinitely differentiable so that its Taylor 
expansion is convergent. If this condition is not met the conver¬ 
gence of the Fourier series is not guaranteed. That is, analyticity 
implies that Dirichlet conditions for convergence are fulfilled. 

In asteroseismogy studies (and often in astrophysics) such 
a convergence of the Fourier series is assumed, i.e. the analyt¬ 
icity of the underlying function is thus also assumed. To date, 
no methods to verify it have been proposed. Here we propose a 
practical method to verify the analyticity condition of the under¬ 
lying function from the observed sampled data, which in conse¬ 
quence provide evidence of the convergence of the Fourier se¬ 
ries. 

In this sense, analyticity refers in this paper to the property 
that a given function can be expanded in Fourier series. To study 
the analyticity of a function we introduce a numerical approxi¬ 
mation for the differentiability of discrete time series: the con- 
nectivitjQ. 


3. Connectivity 

Time-dependent astronomical data are given as a discrete se¬ 
quence of samples from an underlying function, which is a rep¬ 
resentation of a given physical observable. This latter is under- 


^ Avoid confusion with graph theory term connectivity llDiestell200.5l) 
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stood as the physical quantity that satisfies a ll the prop ertied 
derived from the Lebesgue measure definition (lJechlll99% . 

The discrete sampled sequence of photometric variations 
of luminosity and/or radial velocity measurements of pulsating 
stars is commonly called light curve (or RV). Let us consider 
here such a series as a finitely close sequence of data points 
which samples a function composed by a set of infinitely close 
data points. Under these conditions, we want to analyse whether 
such a light curve fully and unambiguously determines the prop¬ 
erties of the sampled function. This is the well-known problem 
of sam pling. It is known (Nyquist-Shannon theorem, IShannonI 
1 19491) that this is possible when the sampled function is band- 
limited to less than half the sampling rate, i.e. the Nyquist fre¬ 
quency. When this happens the complete information of the con¬ 
tinuum (i.e. the sampled function) is contained in the discrete se¬ 
quence. This implies that all the properties of the function can be 
studied, even those concerning the pointwise limits of the func¬ 
tion. Consequently, we are allowed to adequately evaluate the 
differentiability of the sampled function. 

Let us now generalize our discrete sequence of data to sam¬ 
ple a function composed by an stochastic sequence plus a deter¬ 
ministic contribution that can have or not a wave-like structure. 
Such decomposition is guarant eed under th e hypothesis of sta- 
tionarity by the Wold theorem (IWoldlll93^ . The random com¬ 
ponent might complicate the differentiability analysis of a time 
series, however the Kolmogorov continuity the orem allow us to 
fully characterize the properties of the function (lRevuzlll999h in 
a similar way as it is done for pure deterministic signals. This 
requires necessarily the topological separability of the function, 
which is applicable under certain technical constraints related 
with the moments of the data differences, particularly variations 
at most as those given by a power law (see Appendix A for a 
discussion on these constraints). When these restrictions are ful¬ 
filled Kolmogorov theorem states that a stochastic process has a 
continuous representation. 

This is a necessary and sufficient condition for the connec¬ 
tivity analysis that we will define in the next section. 

3.1. Definition 


where At is the sampling rate of the sequence. Note that connec¬ 
tivity is closely related to derivability. In fact, we can define a 
new quantity £)„ as: 




"t" yn+l yn-1 
2Af 


(4) 


which reduces to the typical point derivative (lBurdenll20I ih for 
discrete data when the non-differentiable component of the 
sequence, is zero. Likewise, the connectivity as defined above is 
a function of that is: 


Cn — C(^n) + £n 


(5) 


From this equation, when connectivities are not zero but an inde¬ 
pendent normal stochastic sequenc^^ they can be identified with 
the random term e„ which is a numerical error. In this case the 
non-differentiable component is zero (^„ = 0) and therefore the 
derivatives are still well defined. When connectivities are cor¬ 
related with the signal the non-differentiable component is dif¬ 
ferent from zero (^„ 4= 0) and therefore a derivative cannot be 
defined. In this case, the function is non-analytic and it does not 
satisfy the conditions under which Parseval’s theorem is demon¬ 
strated, meaning that Fourier analysis is not consistent. 

In summary, sampled data are obtained through the sampling 
of functions that needs not be necessarily differentiable a priori. 
Connectivities can be interpreted as jumps in the derivatives giv¬ 
ing rise to discontinuities. This makes connectivities mathemati¬ 
cally self-consistent. They can be considered a s equiva l ent to the 
non-differentiability coefficient introduced bv I Wiene'3 (Il923h in 
a different physical context, and allowing to check the differen¬ 
tiability of the underlying functions. 


3.2. Method 

The forward and backward extrapolations necessary to obtain the 
connectivities have to be calculated numerically. Two conceptu¬ 
ally different numerical approaches were considered. A compar¬ 
ison of both approaches allows us to confirm the results without 
ambiguity. 


Given a data point from a discrete sequence, we can calculate a 
forward and backward extrapolation from the subsets bracketing 
the selected data point. Then, it is possible to check whether both 
extrapolations converge to the same value and coincide with the 
value of the selected data point. 

The above theoretical framework permits the following self- 
consistent definition of the connectivity C„ of a data point y„ of 
a discrete sequence which samples the function F(t) 

Cn -4-4, (1) 

where ef, e^, are the deviations of the forward and backward 
extrapolations from the sample y„ 

4 -yi-yn 

4 ^yn-yn- ( 2 ) 

represented by y{ and y*, respectively. In fact, a discrete approx¬ 
imation of the derivability condition of F at the point t„ can be: 

yj-ynrl ^ yi - yn-l ,,, 

At At ^ 

^ We want to stress here that an observable of a given physical property 
is not necessarily coincident with the function which represents it. 


3.2.1. Cubic splines 

. The Stone-Weierstrass theorem (lRovdenlll988h states that if a 
function is uniformly continuous in a closed interval it can be ap¬ 
proximated arbitrarily well, i.e. as closely as desired, by a poly¬ 
nomial of degree n, n being a natural number. Then, the entire 
function can be represented by a piecewise polynomial func¬ 
tion of finite degree. But piecewise polynomial functions have 
a drawback, the Runge’s phenomenon, showing oscillations at 
the edges of the approximated function. This is equivalent to the 
Gibbs phenomenon in the Fourier approximation. 

Taking all of this into account, we selected piecewise smooth 
polynomials for our modelling. 

Formulated in this general way any particular property can 
be extended to any analytic function. That is, if a continuous 
function can be fitted arbitrarily well with a given cubic spline 
parametric model then the analyticity condition is automatically 
satisfied, otherwise the convergence of the Fourier series is not 
guaranteed for that particular function. In short, the cubic splines 
approach gives a key information about local features of the 
function that generates the observed discrete sequence. 

* This is the theoretically expected distribution of independent sam¬ 
ples. It is usually understood as gaussian white noise. 
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To calculate the coefficients of t he cubic spl ine functions a 
tridiagonal linear system is solved (lBooii[l97^ . The only pa¬ 
rameter affecting the spline interpolation that can be adjusted is 
the number of datapoints of the fitted segments. 

3.2.2. ARMA 

These models dBoxlll9'^ give a parametric representation of 
a time series as a non-closed recursive formula. Autoregressive 
modelling ha s been used be fore for time series analysis in aster- 
oseismology (lKovacslll983h although it has not become popular 
due to the difficulties for the physical interpretation of the fits, in 
contrast to the well-known Fourier frequencies. 

We are interested here in the presence of any non-analytic 
content of the underlying function describing the observed sig¬ 
nals. In contrast to the cubic splines, ARMA models are able to 
represent non-analytic functions, thereby allowing us to identify 
non-analytic signals from the residuals. 

The algorithm used to calcula te ARMA models is sim ilar to 
gap-filling algorithm MI ARM A dPascual-Granadol l20 1 -5h . This 
is based on minim izing rms resid uals and hence, is insensitive to 
the arrow of time dScargIel|1990t) . However, here we extrapolate 
only one point each time and the segments of the modelled data 
are short enough to be able to estimate derivatives without bias. 
The invariance to the arrow of time allowed us to calculate back¬ 
ward extrapolations simply as an inverted forward extrapolation. 

Like cubic splines modelling, the accuracy of the extrapola¬ 
tions with ARMA models depends on the number of data points 
of the segment to be modelled. In addition, ARMA extrapo¬ 
lations depend on two parameters: the orders p and q, which 
are the number of autoregressive coefficients and the number of 
moving average coefficients, respectively. Hereafter, ARMA de¬ 
scriptions of the time series will be denoted by its coefficients as 
ARMA(/7, q). In general, an ARMA model is able to reproduce 
a sinus wave with only two terms. Accordingly, the number of 
terms increases with the frequency content of the signal. The 
ARMA modelling requires of an iterative process involving the 
following steps: 

1. Identification of the order of the starting model using 
Akaike criterion ( Akaikd [19741) and described in detail in 
iPascual-Granadol ( 20I5i) . This affects principally the veloc¬ 
ity of convergence but does not modify the final result. 

2. Calculati on of the para meters using a Steiglitz-McBride al¬ 
gorithm (ISteiglitzill965h given the orders for the AR and MA 
contributions. 

3. Evaluation of the validity of the model based on the good¬ 
ness of the fit in backward and forward extrapolations. The 
algorithm starts again with increased orders until a minimum 
in the residuals is reached. 

4. Connectivity analysis applied to stellar light 
curves 

The test was applied to the light curves of two delta Scuti 
stars observed by CoRoT and Kepler satellites: HD174936 and 
KIC006187665. The parameters required for this modelling are 
specific for each of the studied time series. However, in order 
to avoid misleading interpretation of spurious numerical effects 
we maintain the same parameters when characterising the con¬ 
nectivities for both time series. The 5 Scuti star HD 174936 was 
taken as a reference. 

For the first step of the ARMA approach, the initial tenta¬ 
tive orders were obtained using autoregressive methods (AR). 


Then we iterated until minimal randomly distributed residuals 
were obtained, fixing as optimal orders p -20 and q - Cubic 
splines and ARMA(20,1) connectivities were then calculated for 
every point in a segment of 1000 samples of the time series using 
subsets of 80 data points (40 forward and 40 backward extrap¬ 
olations). We adopted this number of samples per segment as a 
balance between computer time consuming and reliability of the 
ARMA modelling. In any case, we checked that the statistical 
properties of both approximations do not change significantly 
for a higher number of datapoints modelled. 

In order to avoid artifacts we have selected data segments 
with no gaps and homogeneous sampling. Also, studied seg¬ 
ments were statistically normalized allowing comparison of the 
connectivity results obtained using different time series. 

4.1. HD 174936 

We first applied the method to the light curve of the multi 
periodic 6 Scuti star HD 174936 using as analytic solution 
the 422 independent fr equencies found in its periodogram 
(iGarcfa Hernande^l2009h . 

These are distributed in a range of frequencies below 
900juHz. This very large range of excited frequencies is not 
predicted by any non-adiabatic model describing opacity-driven 
pulsators. 

The numerical methods described above were applied to an 
analytic model of this star. The time series of this model was 
calculat ed using the ampli t udes a nd phases of the frequencies 
found in iGarcfa Hernand^ (120091) for the signal component, and 
a noise component based on additive gaussian white noise with 
a level obtained from the same reference. 

For both the observed data (Fig. 1) and the time series of 
the analytic model (Fig. 2), the connectivities calculated using 
cubic splines have higher dispersion than the ARMA ones. Nev¬ 
ertheless, a striking effect is seen in the observed data that can 
not be noticed in the analytic model: spline connectivities are not 
randomly distributed but strongly correlated with the signal, sug¬ 
gesting that these data are not sufficiently well described by the 
analytic model. The same effect is present, although five times 
smaller in amplitude, for the ARMA results. 

The difference between the analytic model and observed data 
has been usually interpreted as correlated noise caused b y the ef¬ 
fects of the turbulence present in the envelopes of stars dChanlinl 
Il997h . In order to check if the convection is at the origin of this 
phenomenon, some additional tests were performed by adding 
non-white n oise to the analyt ic model following the procedure 
described in iKallingeJ (l201f)l) . A power law with two compo¬ 
nents (36 ppm and 19 ppm of amplitudes) was added in the fre¬ 
quency domain to the analytic model (see Appendix B). The re¬ 
sults of this test do not significantly change the previous ones, i.e. 
correlated connectivities were found. We conclude that t he fully 
devel oped turbulence as represented by Harvey models (iHarvevi 
Il985h is not at the origin of this phenomenon. 

4.2. KIC006187665 

In order to discard possible instrumental effects specifically 
linked to a given instrument, a similar test was applied to the 
time series from another star KIC 006187665 observed by the 
Kep ler satellite, and cla ssified as a hybrid Gdor/Dscuti pulsator 
by (IUvtterhoevenll20fl1) . The standard Fourier analysis of the 
light curve supplied by the satellite at a sampling rate of 60 sec¬ 
onds (short cadence regime), yielded 659 significant peaks. 
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Fig. 1. Connectivities of the CoRoT data for the 6 Scuti star HD 174936. a, ARMA connectivities, b, cubic splines connectivities, c, the original 
light curve. (Note the different scaling factors in top and middle panels.) 


The connectivities corresponding to the Kepler target show 
a very similar behaviour to the CoRoT ones (Fig. 3a,b), that is, 
cubic splines connectivities show a more dispersed distribution 
than the ARMA connectivities, and strongly correlated with the 
original Kepler time series. We have calculated the Pearson cor¬ 
relation coefficients for both splines and ARMA connectivities 
of the two time series studied here confirming these results (see 
tabled). 

The non-random structures of the connectivities reflect the 
non-analytic component of the studied function. In order to char¬ 
acterize them, we calculated their amplitude spectrums (Fig. Cl) 
after verifying its applicability by using the test described in this 
paper (see Appendix D). The low frequency range in the case 
of splines connectivities is very similar to that found in the am¬ 
plitude spectrum of the original time series. However a regular 
pattern can be seen along the frequency axis increasing in am¬ 
plitude with frequencies. The ARMA connectivities show some 
pattern at low frequencies whereas an almost flat spectrum, ap¬ 
pearing as a white noise spectrum, is obtained at higher frequen¬ 
cies. The amplitudes of the ARMA connectivities are 1-2 orders 
of magnitudes lower than those calculated with cubic splines. 


5. Discussion and conciusions 

The test introduced here through the connectivities allows us to 
fully characterise the differentiability of the underlying function 
of a time series. We remark that although the underlying function 
can be non-analytic due to the noise itself being non-analytic, 
our test assumes the non-analyticity of this component and eval¬ 
uates the deterministic component through the correlation of the 
connectivities. 


We have shown that at least for the stars studied here the 
functions underlying the observed time series are non-analytic. 
This comes from the finding that connectivities theoretically ex¬ 
pected to be zero or a white noise distribution are strongly corre¬ 
lated with the signal. We have shown this to be an intrinsic prop¬ 
erty of those functions, which implies that their Fourier expan¬ 
sion convergence is not guaranteed. In addition, the fact that pe- 
riodograms actually converge to a real value does not guarantee 
that this value is an asymptotically {N —» oo) unbiased estimator 
of the spectral density of the underlying function. Therefore, the 
value estimated for the spectral density at a given frequency bin 
is not necessarily the true value. 

In order to illustrate the possible relationship between the 
non-analyticity and the unexpected huge number of frequencies 
found in multi periodic stars, we performed a connectivity and 
frequency analysis of simulated dat a from a known function, the 
Weierstrass function (lBailevll2007h . This function is originally 
defined by an infinite sum of harmonic components (see Ap¬ 
pendix E) b ut can be represe nted using eg. IE. 21 with just a few 
components (IWeissteinl2015h . However, it is considered a patho¬ 
logical function since it is continuous everywhere but differen¬ 
tiable nowhere. We show that when performing frequency anal¬ 
ysis of the simulated data the residuals of the prewhitening pro¬ 
cess does not converge to white noise. Notice the large number 
of frequencies obtained as compared with the few components 
used to generate the function (see Eig. lE.3b . This is an example 
of an inconsistency in the application of harmonic analysis that 
can be related to the non-analyticity of the underlying function. 

We are not claiming that Weiertrass function is at the ba¬ 
sis of the inconsistencies found in the time series of the two 6 
Scuti stars observed from space with different instrumentation, 
we are only using its pathological properties to try to "simulate" 
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Fig. 2. Connectivities of an analytic model of the star HD 174936 based on the first 422 frequencies detected using standard Fourier techniques, 
a, ARMA connectivities, b, cubic splines connectivities, c, the original light curve. (Note the different scaling factors in top and middle panels.) 
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Fig. 3. Connectivities of the Kepler data for the hybrid star KIC 006187665 and the original light curve for comparison, a, ARMA connectivities, 
b, cubic splines connectivities, c, Kepler light curve for comparison. (Note the different scaling factors in top and middle panels.) 
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the real cases here studied. Although this peculiar function sug¬ 
gests a Fourier expansion in its standard form (eq. EB Fourier 
frequencies, given by an integer number of the argument, can¬ 
not be defined. This could points towards a reinterpretation of 
the term frequency in the context of the harmonic analysis of 
pulsating stars. The consequences of this new description of fre¬ 
quency are obviously out of the scope of this paper. Although the 
origin of these inconsistencies is unclear, our next step is to try 
to correct numerically the time series in order to overcome the 
analyticity test and therefore obtain a cleaner power spectrum. 

We conclude that periodograms do not provide a mathemat¬ 
ically consistent estimator of the frequency content for stellar 
variability of the objects here studied. Thereby, the concept of 
unambiguous detection of a physical frequency should be re¬ 
vised. This constitutes the first counterexample against the cur¬ 
rent paradigm which considers that any physical process is de¬ 
scribed by a continuous (band-limited) function that is smooth 
and infinitely differentiable. 

This could be on the basis of the unknown nature of many 
phenomena related to the power spectrum of AF pulsating stars: 
the unexpected huge number of frequencies found in multi pe¬ 
riodic stars and the range of these frequencies as mentioned be¬ 
fore, and also the ubiquitous presence of correlated noise in the 
residuals of the fitting of the light curves. 

This paper shows an inconsistency in the application of har¬ 
monic analysis to some pulsating stars observed by space mis¬ 
sions like CoRoT or Kepler. This inconsistency is related with 
the non-analyticity of the underlying function, the origin of 
which is still to be unveiled. 
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Appendix A: Kolmogorov continuity theorem 

As mentioned in Sect.3 this theorem (lRevuzlll999h allow us to 
fully characterize the properties of a stochastic function through 
its continuous representation. The theorem comes after the 
following definition: 

Def.: An extension (?,) of an stochastic process (T,) is a 
process such that for every f > 0, P(Y, = ?,) = 1. As a corollary, 
if the process (Y,) satisfies for all times T > 0 and 0 < s,t < T, 

E(||T,-T,ir)<c|f-sr^", (A.l) 

with a, s, c positive constants, then there exists an extension of 
the process (Y,) that is a continuous process whose paths are al¬ 
most surely continuous. This is known as topological separabil¬ 
ity. 

The conditions above mentioned are quite mild and they al¬ 
low us to consider stochastic processes as continuous paths in¬ 
stead of discrete values as it is the case for deterministic func¬ 
tions. 

It can be demonstrated that these conditions are fulfilled by 
any white noise process and even for a non-gaussian noise hav¬ 


ing a spectral density proportional to /" (with a < 0), i.e. the 
so-called coloured noise. 

As an example, a real-valued Wiener process (so-called 
coloured noise proportional to v“^) has continuous trajectories 
a.sQby this theorem. This means th at it admits a s eparable and 
progressively measurable extension (ICar)asso|l2005h . 

In any case, in this work we refer to the stochastic component 
of a time series with the properties of gaussian white noise. 


Appendix B: The analytic model construction 

The numerical methods described in the main text were applied 
to an analytic model of the 6 Scuti HD 174936. This model was 
built using the amplitudes and phases of the first 422 frequen¬ 
cies fo und using a standard Fourier analysis in iGarcfa Hernand^ 
(l2009h . fitted to an harmonic Fourier-like expansion of the form 

n 

x^ti) = Yj n Vjti + (Pj) +N'^ + A™ (B. 1) 

!=1 

where A,-, vj, and (pj represent respectively the observed ampli¬ 
tudes, frequencies, and phases. The A* -H N'^ terms represent the 
white noise, defined as normalised distribution scaled to the data. 


A* =< > (B.2) 

and the non-white noise, calculated as 

A"w=T^i[Anw] (B.3) 


where A‘™ is defined in the transformed domain as 


A™ 




1/2 


e 


iG 


(B.4) 


where and bk are free parameters. The values of c^t = 4 and 
m - 2 were adopted following Kallinger & Matthews (2010). 
This noise was calculated setting the 0 phases equal to a random 
uniform distribution in the [-tt, n} interval. The light curve so 
constructed is supposed to be originated by the oscillation modes 
of the pulsating star plus gra nulation noise modelled as a fully 
developed turbulence regime (lHarvevll985h and given by eq. B4. 


Appendix C: Pearson correlation coefficients 

The correlations between the connectivities and their corre¬ 
sponding original time series are evident in Fig lC.il even so, 
we calculated their corresponding standard Pearson coefficient 
in order to quantify the correlation between these two series (see 
Table Oil. We also calculated P-values for testing the null hy¬ 
pothesis (i.e. no correlation). P-values were calculated using a 
t-Student distribution. Considering the negligible values ob¬ 
tained for the time series of the selected stars, the null hypothesis 
must be rejected under any computable limit of validity, i.e. the 
two series are strongly correlated. On the other hand, for the case 
of ARMA {Pa), this limit can be put as an standard choice for this 
parameter, e.g. P < 0.05, i.e. only in 5% of the cases the ARMA 
connectivities would be correlated with the original signal. 

Since ARMA models can capture the non-analytic compo¬ 
nent of the signal it is expected such non-correlation of the con¬ 
nectivities. 

^ almost surely in mathematical jargon 
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Fig. C.l. Amplitude spectrum of the connectivities of the Kepler data for the hybrid star KIC 006187665. From top to bottom: panel (a) shows the 
amplitude spectrum of the ARMA connectivities, panel (b) amplitude spectrum of the cubic splines connectivities, panel (c) amplitude spectrum 
of the original light curve for comparison. 


Table C.l. Pearson correlation coefficients ip) and the associated prob¬ 
abilities of no-correlation (P) for ARMA and splines connectivities for 
the time series considered in this work. Subscripts a and s refer to the 
ARMA and splines approaches, respectively. 


Time series 

Pa 

Fa 

Ps 

P. 

CoRoT 7613 

-0.0075 

0.8139 

0.4243 

5.8247e-45 

Analytic model 

0.0072 

0.8199 

8.371e-4 

0.9789 

KIC 006187655 

-0.0257 

0.4176 

0.5553 

5.8039e-82 


Appendix D: Second-order connectivities 

Analyticity, understood as infinite differentiability, guarantees 
the convergence of the Fourier expansion of a function, so al¬ 
lowing to perform Fourier analysis. The test described in this 
work indicates whether the function underlying a time series 
is differentiable or not. Considering the connectivities as a new 
time series, the calculation of the connectivities of the connec¬ 
tivities (from now on second-order connectivities) informs us 
about their differentiability. If, using Pearson cotTelation coeffi¬ 
cients, second-order connectivities are independent, normal, and 
randomly-distributed, then connectivities are differentiable. It is 
thus legitimate to perform a Fourier analysis of them. 

We applied this method presented to the splines connectivi¬ 
ties of the stellar light curve of KIC 006187665 (Fig. Dl, panel 
a). Contrary to the connectivities, second-order connectivities 
calculated using the splines approximation (Fig. Dl, panel b) 
show a normal white noise distribution. These, when calculated 
with splines show that the connectivities are differentiable. The 
ARMA approximation does not provide better fits than the cubic 
splines in this case. 


Appendix E: Non-analytic model 

Frequency and connectivity analysis has been applied to a non- 
analytic model based on a numerical realizatiorQ of the Weier- 
strass function. This is an extreme case of a function exhibiting 
non-analyticity because it is continuous everywhere but differ¬ 
entiable nowhere. 


Appendix E. 1: Definition 

The original definition of the Weierstrass function is: 


W(x) — ^ a’'cos{b"7Tx) 0 < a < 1, 

n=0 


a-b > \ + —n (E.l) 


and b, a positive odd integer number. For rational numbers x - 
—, which is our case, W(x) can be calculated with a finite sum 

q 

such as: 



sin(^) 


(E.2) 


This formula allows to generate an exact numerical realization 
of the Weierstrass function without calculating an infinite sum. 

Eor the non-analytic model that we used we generated three 
realizations of W{x) with different q values and obtained a 30000 
datapoints time series by iterating on p and adding each q com¬ 
ponent. 


Appendix E.2: Connectivity anaiysis 

We apply the numerical techniques introduced before to the time 
series generated as described in previous section. Connectivities 

^ obtained by sampling the continuous function 
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Fig. D.l. Second order connectivities calculated using the cubic splines approximation for Kepler data of the the hybrid star KIC 006187665, 
and the first order connectivities for comparison. From top to bottom: Panel (a) first order spline connectivities, panel (b) second order splines 
connectivities. Note that the correlation present in the first is lost in the last. 


3-component Welerstrass function 



Fig. E.l. A realization of the Weierstrass function with three compo¬ 
nents corresponding to three different q values as indicated in eglOl 


calculated using the splines approximation are clearly correlated 
(Fig lE.21 i and ARMA connectivities are not. These results are 
consistent with the definition and properties of the Weierstrass 
function and demonstrate with the analytic model discussed in 
Appendix B the validity of the connectivity analysis as a test for 
the analyticity of the underlying function of a time series. 

Appendix E.3: Frequency analysis 

A power spectrum of the time series of the non-analyti c model 
and f requency detection was performed using sigspec (iReegenI 
l2007h in order to check the effect of the non-analyticity in a stan¬ 
dard Fourier analysis. A sequence of power spectra was obtained 
during the prewhitening process of different sets of frequencies 
(Fig lE.31 l. The power spectrum at each step was significantly dif¬ 
ferent to a white noise spectrum and the prewhitening process 
continued until the program stopped at 738 frequencies when 
the limit in significance was reached. 


In conclusion, a frequency detection procedure using stan¬ 
dard Fourier analysis is not convergent. This is due to the fact 
that, in this case, the power spectrum calculated with a DFT is 
not a consistent estimator for the original frequency content as a 
consequence of the non-analyticity of the underlying function. 
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Fig. E.2. Connectivities of a realization of a 3-component Weierstrass function. Upper panel: ARMA connectivities. Middle panel: cubic splines 
connectivities. Lower panel: the original light curve. (Note the different scaling factors.) 



Frequency [arbitrary units] 

Fig. E.3. Power spectrum of the time series generated with a realization of a 3-component Weierstrass function (top panel) and the power spectra 
obtained during different steps in the prewhitening process (succesive panels). 
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